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Abstract 

A fundamental non-classical fourth-order partial differential equation to describe 
small amplitude linear oscillations in a rotating compressible fluid, is obtained. The 
dispersion relations for such a fluid, and the different regions of the group and phase 
velocity are analyzed. 
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1 Introduction 



In this paper we study the dynamics of a rotating compressible fluid, taking 
into account such factors as non-homogeneity, and the presence of Coriolis 
forces due to the Earth rotation. The work is devoted principally to the study 
of one of these aspects of the problem, namely the influence of the rotation 
over the wave propagation in the fluid. 

Much work has been devoted to theoretical, mathematical, and physical as- 
pects of waves in compressible and incompressible fluids. The literature of this 
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subject is enormous, primarily because it provides a basis for a wide range of 
applications and techniques, including the solution and modeling of Geophys- 
ical situations, and in the theory of stars. For example the papers [l]-[4], and 
books which discuss the fundamental equations governing waves in such fluids 
[5]. More particularly, study of problems of wave propagation in a rotating 
compressible fluid have been the subject of several earlier publications [6] - 
[7]. Some works have been devoted to the general theory of equations in strat- 
ified and rotating fluids [8] - [10]. Some monographs are devoted to this type 
of problem [11] - [13]. Nevertheless, with few exceptions [14], the systematic 
analtytic study of concrete non-stationary problems of wave propagation in a 
rotating compressible fluid have not been treated. 

The present line of investigation has its origin in the seminal work of S. L. 
Sobolev [15]. Starting from the hydrodynamic equations we obtain in Sec. 2 a 
non-classical fourth-order partial differential equation, which describes small 
amplitude linear oscillations in such a fluid. In Sec. 3 we analyse the disper- 
sion relations for linear waves in the fluid, and discuss the different numerical 
regions of the group and phase velocities of the propagating waves, includ- 
ing their dependence on the direction of the wave vector with respect to the 
rotation axis of the fluid. 

The results obtained are potentially useful for the study of the structure of 
the hydrophysical fluid fields in oceans and in the Earth atmosphere. 



2 Basic equations 



Consider an ideal compressible fluid that fills the whole space. We assume 
that the fluid rotates with constant velocity a/2 around a given axis. The 
fluid motion is refered to a system of Cartesian coordinates x = (xi,x 2 ,x 3 ) 
which rotates with the fluid. The Ox 3 axis is directed along the rotation axis. 

For such conditions, following ideas developed in [11] and [12], the system of 
hydrodynamic equations which describe the fluid is given by the following: 

(1) The Euler equation: 

^ + (v-V)v + -VP-f = 0, (1) 
ot p 

where v = (vi,v 2 , v 3 ) is the vector velocity of the fluid particles, P is the 

— * 

pressure, p is the density of the fluid, and / is the force per unit mass 
acting on the fluid. 
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(2) The continuity equation 

f t+P (V-v)=0. (2) 

(3) The Thermodynamic state equation 

P = P(s,t), (3) 
where s is the entropy of the fluid. 

If we denote (j^) by c 2 , that is to say 

-(f)/ 

then we have 
dP _ 2 dp 



where c 2 is a function of P and p, and c has the meaning of the sound velocity 
in the fluid. Taking into account (5), Eq. (2) can be written 

if + pV-^O (6) 



for iso-entropic motions in the fluid. In (1) / is the Coriolis force due to the 
rotation of the fluid. This means 

/ = -a x v. (7) 



The corresponding Euler equation now becomes 

r\ — * -i 

^ + {v-V)v + -VP + axv = 0. (8) 
at p 



Eq. (8) is nonlinear. Therefore, the exact theory of physical processes in the 
fluid is nonlinear. In this paper, we study a linear approximation, and, to 
this end, we consider small perturbations in the fluid. We therefore linearize 
equation (8). The nonlinear term in (8) is the inertial one, (v ■ V)v. We let 
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vq be the velocity amplitude, 7 the frequency of the wave, and k the wave 
number. Then, we have 

= (v-f)v = kvl = vo_ 

f 7^o V 



where v p h = j/k is the phase velocity of the wave, e is a parameter of nonlin- 
earity. We assume £ < 1, that means that the velocity of the fluid particles 
on the wave is much less than the phase velocity of the wave. Under these 
conditions, we can eliminate the inertial term and obtain 

r\ — * -1 

^ + axv + -VP = 0. (10) 
at p 



This equation, along with (6), comprises the system of equations we shall 
study. We consider small displacements of the pressure and of the density 
with respect to an equilibrium values: 



P = Po + P, P = Po + p', P < Po, p' < Po- 
Under these conditions (10) gives 



dv 1 -±, 

— + a x v + — V(p + p = 0, 

at po + p' 



and (6) gives: 

1 d(p + p) 
c 2 dt 



+ (p + p')V -v = 0. 



By neglecting p' in these expressions, we arrive at the linearized system of 
equations used to describe the small-amplitude motion in the fluid: 

dv 1 -l , . 

— + a x v + — Vp = 0, (11) 
at po 

^|+P»V-ff=0, (12) 



where p is the dynamical displacement of the pressure from the equilibrium 
position p . System (11) - (12) is a 4-equation system of first order, with 
4-unknown variables, p, vi, v 2 , and v 3 , which describes the small- amplitude 
motion (acoustic waves) in the ideal rotating fluid. 
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We will consider what we will call 2D motion in the fluid: by this we mean 
variation of the dynamic pressure and of the velocity for which 



dp 

dxo 



0, 



and 



dv 
dxo 



0. 



From the physical and geometrical points of view, these kinds of motion are 
possible only in domains with infinite cylindrical shape, with the generatrices 
parallel to the Ox 2 axis. In what follows, we suppose that the main value of 
the density is given by p — 1. Then, the system (11)-(12) in the variables v\, 
v 2 , i>3, and p becomes 



dvi dp 
— - av 2 + ^— = 0, 
dt ox i 



dv 2 
~dt 

dv 3 



+ avi = 0, 
dp 



+ 



0, 



dt dx 3 

1 dp dv\ dv 3 
H H - 

c 2 dt dxi dx 3 



0. 



(13) 
(14) 
(15) 
(16) 



From the system (13) - (16) it is straightforward to obtain the following equa- 
tion 



d 2 



1 d 2 u 
^~d^ 



VoM + — U 



a 



,d 2 u 
dx 2 



0, 



(17) 



where V| = d 2 jdx\ + d 2 jdx\ is the 2D Laplacian. Eq. (17) is satisfied by 
the pressure p and also by the components of the velocity v±, v 2 , and v 3 . This 
is a non-classical fourth order partial differential equation, which governs the 
waves in a rotating compressible fuid. 

From (17) we can see that, when a = (which means that the fluid does not 
rotate), the operator L has the form L = ' ! 



-W 2 u 



and therefore the 



solution must coincide with the solution of the 2D wave equation. This allows 
us to conclude that, as a result of the rotation of the fluid, some terms appear in 
the equation that take into account the rotation. The corresponding solutions 
reflect this fact by their dependence on the angular velocity a. Moreover, 
when a — > 0, the solutions of reduce to the well-known solutions for wave 
propagation in a compressible medium. 
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3 Dispersion relations 



In what follows we shall adopt units in which c = 1. Thus 



L[u] 



d 2 \d 2 u 



— V\u + a 2 u — 




(18) 



dt 2 dt 2 



We consider a solution of the form 



U = UqC 




(19) 



where k — (ki, k^) is the wave vector, x = (xi, X3), and 7 is the wave frequency. 
We obtain: 



7 2 (7 2 



a 2 ) - 1 2 k 2 + a 2 k 2 cos 2 = 0, 



(20) 



where 9 is the angle between wave vector k and the Ox 3 axis. Eq. (20) is the 



dispersion relation for plane waves in a rotating compressible fluid. We have: 



We can conclude that there exists a fundamental difference between the prop- 
agation of plane harmonic waves in a rotating fluid and those in a fluid at 
rest. From (21) with a^Owe obtain k — \ that is to say, the well-known 
relation for plane waves in a fluid at rest. We note the dependence of k on the 
value of angle 9 (the direction of the wave vector k). Two special cases are 

— * 

qualitatively different. In the case in which k is parallel to the rotating axis, 
(21) gives k — I7I, which means that in the Ox^ direction plane waves with 
any frequency 7 propagate, including step-like waves (7 = 0). When we have 

— * 

a wave vector k with direction oblique to the rotation axis (0 < 9 < n/2), we 
infer that waves propagate only for certain values of 7: there is a forbidden 
zone for the frequency 7, namely a cos 9 < 7 < a. In the case where k is 
perpendicular to the rotating axis {9 = n/2) we have that only waves with 
frequency 7 > a propagate. 

By differentiating the dispersion relation (20) we obtain the group velocity of 
the waves, thus 




(21) 




(22) 
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Meanwhile, we know that the phase velocity of the waves is 

c = v » = kk = v- (23) 



From (22) we conclude the following: 

(1) For the fluid at rest (a = 0) and, hence, k 2 = j 2 . We obtain 

v g = -k = v ph . (24) 

7 

Thus the group velocity is equal to the phase velocity; there is no 
dispersion, and the energy propagates in the direction k. 

(2) In general, for «^0, the direction of the group velocity v g (that is to say, 
the direction of propagation of the energy [12]) does not coincide with 
the direction of the vector k, and wave dispersion takes place. 

— * 

When the direction of k coincides with the rotating axis Ox s , i.e. when = 0, 
for the group velocity we have 

= 2 72 - h - 7 2 (°' 1 ^r k ) = ^ = ^ (25) 



since k\ — 0, ks — k, and k 2 = 7 2 in this case. Hence, the waves again 
propagate without dispersion. 

If 9 — 7r/2, and under the propagating conditions 7 > a, we have from (21) 
k = V7 2 — ol 2 . Therefore, for the group velocity, we have 



v = (h 0) = -k. (26) 



Meanwhile from (22), 



? ph = -r—£. (27) 

7 Z — a z 



This means that v 9 < f p h and, therefore, we have in this case a normal disper- 
sion of the waves. This is in complete agreement with the fact that, for 9 = n/2 
and k = ^J^ 2 — a 2 , the waves disperse. For | 7 | < a there is a damping of waves 
travelling perpendicular to the Ox 3 axis direction. 

In the general case we can see, using (20), that the expresion 2 7 2 — a 2 — k 2 in 
the denominator of (22) is always a positive number. From (22) it follows that 
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the group-velocity vectors, corresponding to all possible plane waves starting 
from the same begining, are inside the limits of a characteristic cone 



a: 3 | > 



\J a 1 — 7 : 



2 



Xi\. 
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The angle between the group vector and the rotating axis Ox^ is 



6„ r = arctan 



r 




4 Conclusion 

The general behaviour of acoustic waves in a rotating compressible fluid can be 
described in a simple formulation using a fourth order non-classical differential 
equation, which reduces to the classical one when the fluid is at rest (no 
rotation). The analysis of the dispersion relation exhibits the existence of a 
forbidden zone for the frequencies of the propagating wave. We have also seen 
that, in the presence of the non homogeneity of the space in which the waves 
propagate (owing to the rotation of the fluid), the dispersion of the waves 
takes place when the wave vector is not in the same direction as the rotating 
axis. In future work we shall study fundamental solutions of the basic equation 
(17) for several interesting problems involving the diffraction of waves at walls 
within a rotating fluid. 
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